SN 2001em: Not so Fast 



F. K. Schinzel^, G. B. Taylor^'^, C.J. Stockdale'^, J. Granot^, and E. Ramirez-Ruiz^ 

schinzelSmpif r . de 

ABSTRACT 

SN 2001em is a peculiar supernova, originally classified as Type Ib/c. About 
two years after the SN it was detected in the radio, showing a rising radio flux with 
an optically thin spectral slope, and it also displayed a large X-ray luminosity 
(~ 10^^ erg s~^). Thus it was suspected to harbor a decelerating (by then, mildly) 
relativistic jet pointing away from us. About 3 years after its discovery the optical 
spectrum of SN 200 lem showed a broad Ha line, and it was therefore reclassified 
as Type Iln. Here we constrain its proper motion and expansion velocity by 
analyzing four epochs of VLBI observations, extending out to 5.4 years after 
the SN. The supernova is still unresolved 5.4 years after the explosion. For the 
proper motion we obtain (23,000 ± 30,000) km s~^ while our 2cr upper limit on 
the expansion velocity is 6000 km s^^. These limits are somewhat tighter than 
those derived by Bietenholz & Bartel, and confirm their conclusion that late 
time emission from SN 200 lem, a few years after the explosion, is not driven 
by a relativistic jet. VLA observations of the radio flux density, at 8.46 GHz, 
show a decay as starting ~ 2.7 years after the SN. Collectively, the 

observations suggest interaction of the SN ejecta with a very dense circumstellar 
medium, though the implied opacity constraints still present a challenge. 

Subject headings: (stars:) supernovae: individual (SN 200 lem), (stars:) super- 
novae: general 
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Introduction 



On September 15, 2001, iPapenkova &: Li I (120011 ) reported the discovery of an appar- 
ent Supernova (SN) in the galaxy UGC 11794, also mistakenly called NGC 7112 (RA: 
21h42m22.9s DEC: 12d29m54s; z=0.019493 J20 00.0). Throughout t he paper we use Hq 
= 70 km s~^ Mpc~^, and = 1 — = 0.72 flHinshaw et al.ll2008l ). for which the mea- 
sured redshift corresponds to a proper distance of Dp = 83.14 Mpc and an angular distance 
of Da = 81.55 Mpc. Early spect ral line data suggested it was of type I b/c (more likely 
Ic) (IFilippenko fc Chornockl 120011) . It was detected in the radio about 2 years after the 
explosion (jStockdale et al.ll2004l ). initially displaying a rising flux together with an opti- 

( 4.9 — 15 GHz). It al so showed a high 



cally thin spectral slope oc t 



1.9±0.4^-0.36±0.16 



10^^ erg s"^ at 0.5-8 keV flPoolev fc LewinI 120041 ). This made it a 



X-ray luminosity of 

good candidate for harboring a bipolar relativistic jet pointing away from our line of sight 
( IGranot fc Ramirez- Ruia 120041 ) . which has by then decelerated to mildly relativistic veloci- 
ties. Such late time radi o emission is expected in some Type lb / c SNe that may be associated 
with gamma-ray bursts (jPaczyhskil |200 ll : iGranot &: Loebll2003l ) or alternatively produce only 
mildly relativistic bi polar outflows, and could naturally produce a rising radio flux in the 
optically thin regime (IGranot et al.ll2002l ). as the jets interact with the circumstellar medium 
(GSM) and are decelerated. 



Later observations, about 3 years after the explosion ISoderberg et al.l (120041 ). showed 
broad Ha (FWHM 40 A = 1800 km s~^) lines that are not typical of Type Ib/c SN. Observa- 
tions using Very Long Baseline Interferometry (VLBI) were propose d to resolve the source or 
detec t proper motion if it is powered by decelerating relativistic jets ( IGranot &: Ramirez-Ruiz 
20041 ). H owever, despite seyeral attempts with various instruments the source has not been 



resolved JParagi et al.lbood : Isietenholz fc Bartellbood. [2007i IStockdale et al.l 120051 ). An al- 



ternative explanation was sugg ested bv lOhugai fc Ghevalieil J2006h . who developed a model 



in which regular Newtonian ejecta from the supernova explosion collides with a dense, mas- 
sive circumstellar shell (CS). 

In addition to previously published VLBI observations, we add the most recent High 
Sensitivity Array (HSA) observation from February 2007 (about 5.4 years after the explo- 
sion). In order to provide the most accurate measurements and a conclusive picture for the 



proper motion, as well as the expansion y e 



ocity of S N 2001em, we use our observatio n in 



combination with those of IStockdale et al.l (120041 ) and iBietenholz &: Bartell (l2005l . 120071 ). 
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Observations & data reduction 



We observed SN 2001em with the HSA, which consists of NRAO's Very Long Baseline 
Array (VLBA), the phased Very Large Array (VLA), the Robert C. Byrd Green Bank Tele- 
scope (GBT), the Arecibo Radio Telescope (AR), and the Effelsberg Radio Telescope (EB). 
The observation was conducted on February 5th, 2007 with a total time of 10 hours at 8.4 
GHz (bandwidth 32 MHz, total recording bit rate 256 Mbit s"^). VLB! data was correlated 
with the NRAO VLBA processor at Socorro. Analysis was done with NRAO's Astronomical 
Image Processing System (AIPS) and Caltech's Difmap for modeling. In addition we also 
used the VLA (in A configuration) observation on February 18th 2007, which observed the 
continuum of SN 200 lem in L, C and K-Band over a total time of about an hour. The 



three previous VLBA/HSA obse rvations from J uly 2004 (IStockdale et al.l 2004 ). November 
2004 jBietenholz k Bartell 120051 ) and May 2006 jBietenholz k Bartell 120071 ) were re-reduced 
in a consistent fashion. A summary of the observations is presented in Table [1] including 
important parameters of these observations. 

The data was corrected for ionospheric delay, Faraday rotation, and parallactic angle as 
well as for Earth's orientation. We used fringe fitting to calibrate the data for group delay and 
phase rate. The three datasets from previous VLBI observations have been treated equally, 
to ensure comparable results. Our VLBI observation of SN 2001em was phase-referenced 
to JVAS J2145+1115 (1.4° away) and J2139+1423 (2.1° away). The switching intervals 
were ~180 s for J2145+1115, -93 s for SN 2001em, and ~39 s for J2139+1423, which 
was observed in addition about every 5th cycle. We used J2145+1423 to check our phase 
reference (J2145+1115) to get a measure for the quality of our phase referencing as well as 
a position fixed to the International Celestial Reference Frame (ICRF). In addition we used 
J2139+1423 to check astrometry and over results for possible motion within J2145+1115. 



3. Results 

3.1. Position and Proper Motion 

We obtained the results for absolute position, thus for proper motion, by fitting the 
fully calibrated and phase-referenced (u,f) plane data, of the observations from 4 epochs, 
to a circular Gaussian model. This is a good fit for a point source, as we did not see any 
resolved extended emission. However, the first VLBI observation from July 2004 does not 
use J2139+1423 as phase reference check as described above, instead it uses J2139+1316 
which is 1.0° away from our pointing center. Furthermore only the VLBA was used in this 
initial observation. 
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We plot the position including errorbars in Fig. [H We find a circle with a radius of 
0.170 mas that encompasses all the observed positions. This suggests that SN 2001em shows 
little or no proper motion. 

To get an estimate for the error of the proper motion, we plotted the positions of 
J2139+1423 phase referenced to J2145+1115 as shown in Fig. [2] and Fig. [3], which shows 
according t o our data a proper motion of 0.45 mas year~^ in the northwest direction. 
Finer et al.l (120071 ) give a relative motion of the 4 components between 0.108 and 0.593 
mas year~^ separated up to 3.5 mas from the core. This can be explained by changes in the 



exten ded emission of J2 139 +1423, which shows an extremely curved jet (ISavolainen et al. 



20061 ). in combination with any proper motion of the phase chec k calibrator. The int rinsic 
motion of J2139+1423 has been determined to be < 40 /xas year~^ (IFeissel- Vernierll2003l ) . We 
also assume our reference source J2145+1115 to be stationary. Using a linear relationship 
between the stationary phase center and the apparently moving J2139+1423, that are 3.49° 
apart from each other, we find an error for the proper motion of SN 2001em of ±0.184 mas 
year~^. 

This high error is not ve ry satisfactory and yielded unreasonably small values. As 
shown by iPradel et al.l (120061 ) . the linear approach is not reliable. Thus we had to find a 
better way to obtain a reasonable err or for the posi t ions. Using the discussion of astrometric 
accuracies of VLBA observations by iPradel et al.l (120061 ) and the value of the weighted 
least squares fit, we were able to reduce the position errors to 0.013 mas and 0.051 mas in 
R.A. and Dec. respectively. Using a weighted least squares fit for the motion in R.A. and 
Dec. separately (Fig. H]), we get an error of 0.037 mas year~^ in R.A. and 0.071 mas year~^ 
in Dec. 

We obtained a proper motion of (0.030 ± 0.037) mas year~^ in R.A. and (—0.053 ± 
0.071) mas year~^ in Dec, which gives a total proper motion of 

/i= (0.061 ±0.080) mas year-^ (1) 

The proper distance to SN 2001em is 83.14 Mpc, then the proper motion gives us a projected 
velocity on the sky of 

V = (24, 000 ± 32, 000) km s'^ (2) 



Bietenholz fc Bartell (120071 ) report a proper motion, based on two epochs, of (0.089 ± 
0.093) mas year~^ or (33,000 ± 34,000) km s~^, respectively, where for the cosmological 
parameters we use the latter becomes (35,000 ± 37,000) km s~^. 
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3.2. Size Limits of the Radio Emitting Region 

Again we used the fuUy cahbrated and phase-referenced to J2145+1115/J2139+1423 
{u,v) plane data, to image the SN using Difmap. Fig. [5] shows the VLBI image of SN 
2001em on 2007, February 5. The source is not resolved. The resolution of our observation 
was (2.02 X 0.858) mas at -7.17°. 

A precise determination of the upper limit for the size of SN 2001em was obtained using 
Monte-Carlo simulation. We again used model fitting in the {u,v) plane with a circular 
Gaussian model. We created 100 {u,v) datasets for a range of sizes for each epoch based 
on flux density and corresponding image RMS noise using AIPS task UVMOD. This creates 
simulated datasets with similar properties to the actual data. We found the minimum for 
each simulated dataset and examined the Gaussian distribution of modelflt sizes depending 
on the simulated source diameter deflned in UVMOD. This has been done for all 4 epochs 
and a range of simulated source sizes. The combined results of the upper limits for size 
(FWHM diameter of the circular Gaussian) with a confldence of 95% are listed in Table [2J 
This means that 95% of the Gaussian distribution is below the upper size limit thus giving 
us a 2a upper limit. 

In Fig. O we show a plot of the upper limit radii vs. time and the linear flt to obtain 
an upper limit for the expansion velocity. Assuming this linear expansion behavior since its 
discovery on 2001.704, we flnd a 2a upper limit for the expansion velocity of 0.015 mas year~^, 
which corresponds to 

Vcxp < 6,000kms~\ (3) 



Bietenholz fc Bartell (120071 ) report a 2a upper limit for the average expansion velocity since 
the explosion of 25,800 km s~^, at 4.6 years after the SN. For that same epoch, when 
reanalyzing the sarne data we obtain a 2a upper limit of 8,100 km s~^. In contrast to 



Bietenholz &: Bartell (120071 ) we used a circular Gaussian model for model fltting instead of 
an elliptical one which already reduces the size limit by a factor of about 2. In addition, we 
flt directly to the visibility data and used the Monte Carlo simulation method which in the 
case of an unresolved source improves results further compared with image plane fltting the 
elliptical Gaussian using the AIPS task JMFIT. 



3.3. Radio Lightcurve and Spectrum 

Many multifrequency observations of SN 2001em have been conducted with the VLA 
since 2003. We present the published values and the observations from February 2007. In 
Fig. [7] we show the 8.46 GHz detections plotted over time since the supernova explosion. 



- 6 - 



From the peak brightness the 8.46 GHz peak luminosity is (1.54 ± 0.10) xlO^^ erg s ^ Hz ^. 
We see a power law decay in the flux densities since its peak around 2004 (or ~ 2.7 yr since 
the SN) with an index of a = 1.23 ± 0.40, i.e. oc t"". 

Fig. [8] shows the spectrum at the time of our observations (t = 2007.094,2007.132) 
where our measurements are at 1.46 GHz, 4.9 GHz, 8.4 GHz and 22 GHz, the values are 
listed in Table [31 The 22 GHz point is a Sex detection. The spectral index [F^ oc at the 
cut-off between 8.4 GHz and 22 GHz is /^drop = 0.38 ± 0.18 and between 4.9 and 8.4 GHz it 
is /?platcau = 0.077 ±0.011. 



4. Discussion 

The lack of any apparent proper motion is consistent with other VLBI measurements 



liar), 1993J (lib), 2001gd (lib), and 2004et (IIP) ( 


Perez- Torres et al. 


2C 


02I;] 


Manchester et al. 


2002; 


Marcaide et al. 


2008; 


Perez- Torres et al. 


2OO5I; 


Marti- Vidal et al. 


200' 


r). With a fourth 



epoch, we have further constrained the value reported to (24,000±32,000) km s~^. While 
such constraints on the proper motion may sometimes indicate a symmetric blast-wave cen- 
tered on the supernova position, in our case the 2cr upper limit on the average proper velocity 
(87,000 km s^^ ~ 0.29c) is much larger than that on the average expansion velocity (6000 
km s~^), and therefore no proper motion is expected to be detected even for a reasonably 
asymmetric supernova (i.e. given the limit on the expansion velocity, we would expect to 
detect proper motion only for an extremely asymmetric explosion, which has not yet been 
observed in ordinary supernovae). 

The cons traint on the proper motion d oes, however, rule out a relativistic jet pointing 
away from us (iGranot fc Ramirez-Ruial2004j ) as the source of the radio (and X-ray) emission. 
In such a scenario, the current limit on the size of the emitting region allows only for a very 
compact emission region, and rules out a double sided jet where both sides are visible (even 
if the emission from each side is from a compact region). Moreover, the limit on the proper 
motion also rules out the possibility that only one side of the jet is visible, since for a 
compact relativistic emission region moving at an angle 9 relative to our line of sight, the 
observed limit on the apparent velocity I3ap = /3sin^^/(l — jScosO) < 0.29 implies, e.g., 
9 < 1.86° = 0.032 rad for j3 > 0.9. Such a small angle relative to our line of sight would not 
produce an optically thin rising flux ~ 2 years after the SN (since by that time the Lorentz 
factor is no more than a few and our line of sight would be well within the beaming cone of 
the radio emission), which was the original motivation for suggesting such an explanation. 
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Comparing results for SN 200 lem to SN/GRB pairs 1998bw, 20031w and 2006aj, peak 
luminosities at first suggest a possible GRB event of SN 2001eni , but its peak late in time does 
not match the pattern of SN/GRB pairs (IKaneko et al.l 120071 ) . In additio n GRB SN show 
much larger expansion velocities, >100,000 km s~^, than what we obtained (jPihlstrom et al. 
20071 ). Moreover, our results rule out even a mildly relativistic outfiow with an energy 



comparable to that of the SN ejecta (which is > 10^^ erg), which would not be capable of 
producing a GRB, and may potentially be much more common in core collapse S Ne compared 
to highly relativistic outfiows of similar energy (IGranot fc Ramirez-Ruig|2004j ). 



Our la limiting value for the velocity of the radio-emitting region of a blastwave of 
~ 3,000 km s~^ is unusually low (see Table Hj) and is close to the lowest measured velocity 
of the radio emittin g region of any SN , repor ted for SN 1987A by [Manchester et al.l (120021 ) 



of ~ 3,500 km s'^ iMarti-Vidal et all (120071 ) report a minimum, measured expansion ve- 



locity of (15,700 ± 2,000 ) km s ^ for the radio emitting region of the type IIP SN 2004et. 
Perez- Torres et al.l (120051 ) reported an upper limit for the expansion veloc i ty for the radio 



emitting region of the type lib SN 2001gd of 23,000 km s~^. iMarcaide et al.l (120081 ) report an 
evolving velocity for the radio-emitting region with the outer velocity of the radio emitting 
shell to be ~ 10,000 km s~^ and the inner velocity to be ~ 7,000 km s~^ occurring between 
days 1063 and 1399 since explosion. The onl y other type Iln SN with a measured veloc- 
ity of its radio emitting region was SN 1986J, iPerez- Torres et al.l (120021 ) report an average 
expansion velocity that dropped from 7,500 km s~^ in 1988 to 6,300 km s~^ in 1999. 

What is of particular interest is that S N 2001em, origina l ly cla ssified as s type Ic SN, 
has such a low average expansion velocity. ISoderberg et al.l (120031 ) estimate that the ra- 
dio emitting region's expansion velocity if similar to the type Ic SN 2003L, should exceed 
16 ,000 km s~^ and likely be as hig h as 30,000 km s^^. Assuming the early optical classification 
of iFilippenko fc Chornockl (12001! ) to be valid, we might expect the blastwave of SN 2001em 
to have been resolved near the limit we establish had it undergone a period of rapid expan- 
sion prior to hitting a very dense GSM. Preliminary models of VLA observations suggest 
an ave rage mass-loss of ^ 2.5 x 10~^ M^;^ year~^ , assuming win d -drive n, clumpy/filamentary 



GSM ( IWeiler et al.ll2002l : IStockdale et al.l 120071 ). IWeiler et al.l ((20021) indicates that such a 
radio-derived mass-loss rate is slightly higher than values determined for other type Iln SNe, 
1.14x10"^ Mq year-^ for SN 1988Z and 4.28 xlO'^ M© year'^ for SN 1986J. However, de- 
pending on when this shell of material was ejected, the actual mass-loss event was likely 
much higher than this average value. Due to a lack of any radio observations within the 
first two years following the explosion, we are unable to place any reasonable upper-bound 
on this mass-loss rate. 



We obtained a 2a upper size limit of the radio emitting region of 0.98 x 10^^ cm after 5.4 
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years, which is marginally consistent with the model of IChugai &: Chevalierl (120061 ). in which 
the ejecta collides with a dense circumstellar shell at a radius of ~ 6 x 10^^ cm after ~ 2.6 yr 
(note that a thin emitting shell, that is more relevant for their model, gives an upper limit on 
the radius and average expansion velocity that is smaller by ~ 10% compared to the circular 
Gaussian that we have used). Our results support the picture that at a very early epoch, 
this supernova has evolved from a type Ib/c SN with little circumstellar interaction to a 
type Iln with strong circumstellar interaction. The c lassification as type I In SN was in part 
due to the late-time observation of a narrow Ha line (ISoderberg et al.ll2004r). This d etection 
indicates the existence of cool shocked circumstellar gas fIChugai &: Chevalierl |2006| ) . 

The large X-ray luminosity of Lx ~ 10^^ erg s~^ at t 950 days implies a radiated 
energy in the X-ray range (0.5-8 keV) of the order of tLx{t) ~ 10^^ erg (assuming the 
typical time for variation in Lx is At ~ t). This sets a strict lower limit on the energy 
that was dissipated and converted into internal energy by that time, of -E'dis(i':) ^ 10^^ erg. 
However, the true value of -Edis(^) is most likely significantly higher than this lower limit, 
due to the combination of various inefficiencies in the conversion of the dissipated energy 
into radiation in the X-ray range within the dynamical time. First, the total radiated energy 
is probably somewhat larger than the measured value, since u Fy is still rising within the 



observed energy range {lyFi, oc z/0-9±o-35j iPooley &: LewinI 120041 ) and should therefore peak 
at higher energies. Second, not all of the dissipated energy is given to the electrons, and 
most of it can go into protons (or other ions; both the thermal population or cosmic rays) or 
magnetic fields, and thus would not be radiated away. Finally, even the energy that does go 
into electrons is not always radiated away efficiently on the dynamical time. Given all the 
above, we consider -Edis(i) ^ 10^^ erg to be a more realistic estimate of the energy that has 
been dissipated by time t (although a somewhat lower value might still be possible under 
some circumstances). This would suggest a deceleration time of tdec ^ 10^ days, i.e. that 
by that time the ejecta had interacted with a s urrounding mass comparable to its own mass 
(iGranot fc Ramirez-Ruijbooi I Chevalierl 120071 1. 

Moreover, the radio light curve also suggests a deceleration time of tdec ~ 10^ days, while 
our limits on the source size imply that the typical initial velocity Vq of the mass that has 
decelerated by tdec is ^ 8, 000 km s~^. For an external density p = Ar'"^ this implies = 
A/(5xl0" gcm-i) > 7.3x103(^/10^^ erg) (tdec/10Mays)-Hwo/8xl03 km s-^)-=^ that corre- 
sponds to a mass-loss rate of M > 0.073(t;^/103 km s-^){E/10^^ erg)(tdec/10Mays)-Ht;o/8 x 
10^ km s^^)^'^ Mq yr^^ where is the stellar wind velocity of the progenitor star, as- 
suming a quasi steady wind during the relevant time scale. Such a high mass loss rate 
may suggest an episodic mass ejection event, rather than a steady wind, that would re- 
sult in a dense GSM shell. The required mass of such a GSM shell is at least compa- 

> Mo 



rable to that of the SN ejecta that was decelerated until tdec, Mshdi 



2E/vl 
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1.57(-E'/10^^ erg)(i;o/8 x 10^ km s~"^)~^ Mq. Now, let us consider the opacity of this part 
of the original SN ejecta, of mass Mq and initial velocity fo, at earlier times. Assum- 
ing one free electron per proton mass, its Thompson optical depth at t < tdec would be 
tt = aTE/{27rmpV^t^) = 23(^/10^1 erg) (wo/8 x 10^ km s-i)-^(t/30 days)-^, while additional 
opacity from other parts of the ejecta could only add to the total optical depth. This 
opacity is somewhat large for the time near the peak of the supernova optical light curve 
(if: ~ 30 days), and may become a more severe problem if vq turns out to be lower than the 
current upper limit on it. 



5. Summary & future observations 

It is clear that SN 2001em exhibits radio properties consistent with a type II SN with 
a significant CSM/blastwave interaction. Future VLBI monitoring is recommended. As the 
blastwave expands we will be able to further constrain or eventually determine the average 
velocity of the blastwave. SN 200 lem has faded considerably at 8 GHz over the last few 
years but is still readily detectable with a estimated flux density at the end of 2008 of 0.55 
mJy. With HSA recording rates expected to increase to 4096 Mbps by 2011, the thermal 
noise level of the HSA at 8 GHz should decrease to ~1.40 yuJy beam~^ at an integration 
time of 8 hours, thus keeping pace with the fading of SN 2001em which is expected to follow 
a power law decay. At the end of 2008 an HSA observation (512 Mbit s~^) would have a 
signal to noise ratio (SNR) of 139 and by 2011 with 4096 Mbit s'^ bandwidth a SNR of 288 
based on la using natural weighting. 
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Table 1. Observational Summary 



Date 


Time after SN 


Frequency 


Integration 
Time 


Bandwidth Polarization 


Instrument 




(days) 


(GHz) 


(minutes) 


(MHz) 





1021 8.410 178 32 2 VLBA 

1165 8.410 400 32 2 VLBA+AR+EB+GBT+Y27 

1686 8.406 423 32 2 VLBA+AR+EB+GBT+Y27 

1969 8.410 300 32 2 VLBA+AR+EB+GBT+Y27 

1983 1.465 12 50 2 VLA 

4.885 9 50 2 VLA 

22.485 20 50 2 VLA 



Notes.- AR = 305 m Arecibo telescope. EB = 100 m Effelsberg telescope. GBT = 105 m GBT. Y27 = phased VLA. 



Table 2. Results for the 2a upper size limit of the FWHM diameter of a circular Gaussian 
point source, obtained from Monte Carlo simulation. The position offsets are given using 
the reference position 21h42m23.619s R.A. and +12d29m50.299s Dec. 



Epoch 
(years) 


Flux Density 
(mjy) 


Oflfset R.A. Oflfset Dec. RMS Sizehmit 
(mas) (mas) (/iJy/beam) (mas) 


2004.496 
2004.893 
2006.401 
2007.094 


1.41 
0.60 
0.52 
0.59 


-0.089 0.15 59 0.19 
0.028 -0.15 36 0.23 
0.083 0.061 21 0.19 
-0.019 -0.057 21 0.16 


Table 3. 


Results for flux densities and corresponding noise of SN 2001em on February 18, 






2007. 


Frequency 
(GHz) 


Flux Density 
(mJy) 


Noise 
(mJy) 


1.465 
4.885 
22.485 


0.704 
0.812 
0.345 


0.264 
0.070 
0.056 



2004 Jul 01 . 
2004 Nov 22 

2006 May 27 

2007 Feb 04 . 
2007 Feb 18 . 



- 13 - 



Table 4. Expansion Velocities for Type II Supernovae. 



Supernova 


Type 


Velocity 
(10^ km s-i) 


Reference 


1987A 


Ilpec 


~ 3.5 


Manchester et al. ('2002') 


2004et 


IIP 


> 15.7 ±2.0 


Marti'- Vidal et al. (2007) 


2001gd 


lib 


23.0 


Perez-Torres et al. (2005') 


1993J (outer shell) 


lib 


- 10.0 


Marcaide et al. (2008) 


1993J (inner shell) 


lib 


- 7.0 


Marcaide et al. (2008) 


1986J (measured 1988) 


Iln 


7.5 


Perez-Torres et al. (2002) 


1986J (measured 1999) 


Iln 


6.3 


Perez-Torres et al. (2002) 


2001em (measured 2006) 


Iln 


< 8.1 


this manuscript 


2001em (measured 2007) 


Iln 


< 6.0 


this manuscript 
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Fig. 1.— Positions of SN 2001em from 2004 to 2007, centered at R.A. 21h42m23.619s Dec. 
12d29m50.299s (J2000). Circle centered at the origin has a radius of 0.170 mas. The triangle 
on the middle right shows the projected velocity vectors in R.A. and Dec. respectively, as 
well as the total projected velocity for a period of 2.201 years (3 epochs of HSA observations). 
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Fig. 2. — Positions of the phase check cahbrator, J2139+1423, from 2004 to 2007. 8 epochs 
of VLBA observations show an apparen t relative motion of its components between 0.108 
and 0.593 mas year~^ (jPiner et al.l 120071 ). 
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Fig. 3.— VLBA images of J2139+1423 (PKS 2136+141) at 15 GHz. The figure includes 
observations from the VLBA 2 cm Survey, from the MO J A VE Survey and a 15 GH z image 
from the multi-frequency data set observed in 2001.36 by ISavolainen et aLl (120061 ). In all 
images, a uniformly weighted (M,u)-grid is used. The Gaussian components fitted to the 
visibility data are shown as ellipses overlaid on each ir nage. The size and orie ntation of the 
beam is shown in the lower left corner of each image (ISavolainen et al.ll2006l ). The images 
show clearly the existence of an extremely curved jet with an inner position angle ranging 
from -30° to -90°. 
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Fig. 4. — Weighted least square fits for position vs. time in R.A. and Dec. This gives a proper 
motion for SN 2001em of (0.030 ± 0.037) mas year'^ in R.A. and (-0.053 ± 0.071) mas year"^ 
in Dec. 
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Fig. 5. — VLBI image of SN 2001em from February 4th, 2007. Peak brightness is 586 //Jy 
bm~^, the rms background noise was 21 //Jy beam~^. The contours are drawn at levels 11, 
25, 39, 53, 67, 81, 95 and times 1% of the peak brightness. The lowest contour is at 3.1a. 
The grey scale flux range is 550 /iJy beam~^. The beamsize is shown in the lower left corner. 
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Fig. 6. — 2(7 upper limits for the radius of SN 2001em plotted over time. The linear curve rep- 
resents the 2(7 upper limit for the expansion velocity of 6000 km s~^ plotted from 2001.7041. 
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Fig. 7. — Plot of the 8.46 GHz lightcurve us i ng our data from February 2 007 and previously 
published values fiBietenholz fc Bartell[2007l : Istockdale et alibooi boosh . The plot shows a 
power law decay in the lightcurve since the peak around 2003. The slope of the power law 

fit is y OC t-l-23±0.40_ 
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Fig. 8. — The observed radio spectrum using the 8 GHz flux density from the phased VLA 
as of February 4 2007 and 1.2 GHz, 5 GHz and 22 GHz observed on February 18th, 2007. 



